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We present a method for parametrizing heavy meson semileptonic form factors using 
dispersion relations, and from it produce a two-parameter description of the B — > B elastic 
form factor. We use heavy quark symmetry to relate this function to B — > D*W form 
factors, and extract \ V c b\ = 0.037^q'qq2 from experimental data with a least squares fit. Our 
method eliminates model-dependent uncertainties inherent in choosing a parametrization 
^ ■ for the extrapolation of the differential decay rate to threshold. The method also allows a 
description of B — > Dlv form factors accurate to 1% in terms of two parameters. 
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1. Introduction 



A nonperturbative, model-independent description of QCD form factors is a desirable 
ingredient for the extraction of Cabibbo-Kobayashi-Maskawa parameters from exclusive 
meson decays. Progress towards this goal has been realized by the development of heavy 
quark symmetry]]]], which relates and normalizes the B — > D*W and B — > Dlv form factors 
in the context of the -A? expansion, where M is the heavy quark mass. This normalization 
has been used|2|-[]]] to extract the value of the CKM parameter \V c b\ by extrapolating the 
measured form factor to zero recoil, where the normalization is predicted. 

This form factor extrapolation, necessary because the rate vanishes at zero recoil, 
introduces an uncertainty in the value of \V c b\ due to the choice of parametrization. Esti- 
mates of this uncertainty obtained by varying parametrizations suffer the same ambiguity. 
This ambiguity could be eliminated if one had a nonperturbative, model-independent char- 
acterization of the form factor in terms of a small number of parameters. 

In this paper we use dispersion relations to derive such a characterization and apply it 
towards the extraction of \V c b\. The characterization uses two parameters that describe the 
form factor over the entire physical range to 1% accuracy. For computational convenience 
we use heavy quark symmetry in our characterization, but other than the normalization 
at threshold, this is an inessential ingredient which may be discarded at the cost of some 
extra algebra. 

In Sec. 2 we describe a well-known method|J for using QCD dispersion relations 
and analyticity to place constraints on hadronic form factors. We then derive a basis for 
functions that obey the constraints imposed on the B — > B elastic form factor F, and show 
that to 1% accuracy, only two terms in the basis function expansion need be kept. In Sec. 
3 we use heavy quark symmetry to relate F to the Isgur-Wise function, which describes 
the form factors for B — > D*W in the infinite quark mass limit. We make a least squares 
fit to CLEO@, ARGUS ||, and ALEPH[| data using \V cb \ and our two basis function 
parameters as variables, and present our results. Reliability of the method is discussed in 
Sec. 4, implications for \V u b\ are discussed in Sec. 5, and concluding remarks are presented 
in the final section. 
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2. The Analyticity Constraints 

Consider the form factor F, defined by 



(B(p')\V^\B(p)) = F(q 2 )(p + p')^ (2.1) 

where = 67^6, and q 2 = (p — p') 2 is the momentum transfer squared. A dispersion 
relation for the two point function (0|TV^T^|0) connects its perturbative evaluation with a 
sum over positive-definite terms. This sum includes a contribution ~ | (0| V^\BB) | 2 which, 
by crossing, is given by the analytic continuation of F. This procedure leads to a bound 0- 
| ona weighted integral of \F\ 2 over q 2 > AM 2 . 

A key ingredient in this approach is the transformation that maps the complex q 2 
plane onto the unit disc \z\ < 1: 



AM Z B 1 - z 



In terms of the angular variable e lB = z, the once-subtracted QCD dispersion relation may 
be written as 2 ^ 

de\^ )F{e^)\ 2 <- , (2.3) 
2tt J 7T 

where the weighing function <p(z) contains both the Jacobian of the variable transformation 

and the essential physics of the perturbative QCD calculation^: 



^ )= iW^ (i+z)Vr3 " (2 - 4) 

Here rif is the number of light flavors for which SU(rif) flavor symmetry is valid; we take 
rif = 2. Perturbative corrections to the dispersion relation are incorporated in p, which 
has been computed [|,[T(J to 0(a s ), p = l + 0.73a s (m b ) w 1.20. 
In terms of the inner product defined by 

(f,g) = ^- Fd0r(9)g(9), (2.5) 
^ Jo 

and the function 

= <j>{e w )F{e w ) , (2.6) 

the dispersion relation Eq. ( |2.3j ) reads 

(*,*)<-. (2.7) 

7T 



Physically, poles of F inside the unit disc originate from resonances below threshold and 
cannot be ignored 0; for the B system these are the resonances Ti 5 2,3- A simple but 
effective trick [§| eliminates the poles with no reference to the size of their residues but 
rather only their positions (i.e., masses). Define the function 

p/ x _ (z- Zl )(z- z 2 )(z-z 3 ) , , 

[Z) ~ (l-z lZ )(l-z 2 z)(l-zszy {2 -*> 

where the Z{ correspond to the values q 2 = My. , and the constant 

T = -P(O)0(O)F(O). (2.9) 

By 6-number conservation, F(0) = 1. Since P{z) has modulus one on the unit circle, the 
function 

f (z) = ±[<f>(z)P(z)F(z)+T\ (2.10) 
is analytic on the unit disk and obeys 

(/o,/o)< --\T\ 2 = I. (2.11) 

7T 

Any function g(z) that is analytic in the unit disc and obeys Eq. ( |2.11| ) may be expanded 
in an orthonormal basis 

oo 

g(z) = Y,anZ n , (2.12) 

n=0 

where a n = (z n , g). The QCD form factor may therefore be written as 
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F(z) 



P(z)<P(z) 
with 



oc 



J2*nZ n+1 ~T 



n=0 



(2.13) 



VM 2 <i-|T| 2 . (2.14) 
z — < n 

n=0 

In the next section we will use heavy quark symmetries to relate F to the form factors 
for B — > D*W, where the physical kinematic range is < z < 0.056. Thus, retaining only 
ao and a\ induces a maximum relative error of vT(0.056) 3 /P(0)^(0) ~ 0.009. 



3. Extraction of \V c b\ 

3.1. Heavy Quark Symmetry Relations 

In the infinite b and c quark mass limit all the form factors for B — > Dlv and B — > Z}*Zi> 
are given by one universal "Isgur-Wise" function. This allows us to apply the constraint 
on F to the particular combination of form factors actually measured, rather than deriving 
constraints for each form factor separately. The Isgur-Wise function is related to the form 
factor F by F(u) = 775^(0;), where u = v-v', and the short-distance matching correction r)B 
is unity at threshold by the Ademollo-Gatto theorem. To compare with data, we need the 
analogous short-distance matching and running correct ion |IT| for B — ► D* form factors. 
For example, for the vector current form factor g one may write g = j^F. This relation 
generally holds to order 1/M, but at threshold it holds to order 1/M 2 |12|]. We treat rj^/rjE 
as approximately constant and equal to 0.985. The errors from this approximation should 
be no larger than the neglected 1/M corrections. 



3.2. Maximum Likelihood Fit 



Once the essential physics of QCD is incorporated into the calculation via Eqs. ( |2.13| ) 
and ( p.!4| ) , the maximum likelihood fit is simply an ordinary chi-squared minimization with 
parameters \V c h\ and the basis coefficients {a n }. As mentioned previously, the smallness of 
z in the physical range allows us to ignore all coefficients except ao and a±. We normalize 
input data to a B lifetime of tb = 1.61 ps. 
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Table 1 Fit values for \V c b\, do, and a± from the various experiments. 



Table 1 shows the central values and 68% confidence levels for \V c b\, Go, and a± from 
the various experiments. The saturation of its QCD bound by a± is not significant because 
its variance is large, which arises because the contribution of a± is suppressed by an extra 
power of z. 
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Figure 1 x 2 per degree of freedom resulting from a least squares fit in \ V c b\, and ao and 
ai from Eq. ($.13Q to CLEO, ALEPH, and ARGUS data, plotted against \V cb \. 



Plots of x 2 P er degree of freedom versus \V c b\ for each of the experiments are shown 
in Fig. 1. The minimum x 2 is consistently low, remarkable agreement for a first-principles 
parametrization. 

Figure 2 shows the product of the best fit form factors with \V c b\, superimposed with 
experimental data. At 90% confidence level, ao and a± are consistent with zero, suggesting 
the dispersion relation may be saturated entirely by higher states. 

The errors on \V c b\ in Table 1 are statistical only; the treatment of systematic errors 
depends both on our parametrization and a detailed understanding of the experiment. 
The error implicit in the variation over choices of parametrization, however, is absent. For 
the ARGUS experiment, varying over four possible parametrizations induced a spread of 
0.012 in \Vcb\i and was the major impediment in using heavy quark symmetry to obtain 
a model-independent extraction. Even the experiment with highest statistics, CLEO, 



remains sensitive to the choice of extraction. For a linear fit, they find \ V c b \ • 10 3 = 35.1 



+ 1.9 

-1.9; 



while for a quadratic fit, they find \V c b\ • 10 3 = 35.3^3 3. Presumably a higher-order fit 
would yield even larger variances. Fortunately, the basis function approach does not yield 
statistical errors indicative of such a higher-order fit. 
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Figure 2 Best fit values for the product of \ V c b\ with the Isgur-Wise function for CLEO 
(solid line), ARGUS (dot-dashed line), and ALEPH (dashed line) data. The data for 
each experiment, adjusted for the B lifetime and zero-recoil normalization used in the 
text, is superimposed. 



4. Discussion of the Method 



4-1. Reliability of the Parametrization 



The parametrization Eq. ( |2.13|) was used to make a three-parameter (ao, ai and \V c i,\) 
maximum likelihood fit to the data. We can list four types of correction to this equation. 

First, perturbative 0(a^(mt,)) corrections to the dispersion relation ( |2.11| ) arise as 
^-independent renormalizations of the function (p, and may be calculated systematically 
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through higher loop diagrams. The extraction of \V c b\ is rather insensitive to such correc- 
tions: Altering by hand the perturbative computation by ±10% changes the central value 
of \V c b\ by less than ±0.3%. 

Second, non-perturbative corrections to the dispersion relation may be analyzed via 
an operator product expansion. The first correction takes the form of a gluon condensate, 
estimated to be^g Q ff^ } < G^G^ >^ 10" 5 , which is completely negligible. 

Third, the truncation of Eq. 13| ) at finite n introduces an error proportional to 
p(o)(j}(o) Si^n+i a i zX+1 an d additionally suppressed by the coefficients a n themselves, 
whose sum is bounded: Yl™=o \ a n\ 2 < I ~ 0.31. Again, for the case at hand, the first 
two terms are sufficient to describe any form factor allowed by the dispersion relations to 
within 1%. 

Finally, the application of the B —>■ B dispersion relation to B — > D*W decays relies 
on heavy quark symmetry, which is computationally convenient but unnecessary. The 
weighing function <p(z) appropriate to B — > DW may be readily deduced from an analogous 
computation for B — > nW[15\. Eq. (|2.13| ) then holds with the function P(z) altered to 



reflect B c poles below threshold, the positions of which have been calculated in the context 
of a nonrelativistic potential model[|TB[. It should also be possible to derive analogous 



dispersion relations for each of the B — * D*W form factors, again using no assumptions 
about heavy quark symmetry. This should allow the extraction of \V c b\ from both B — ► 
D*W and B — > DW decays, given only the normalization of the form factors at zero recoil. 

Alternatively, one may want to use heavy quark symmetry to relate the B — ► D*W 
and B — * DW form factors and construct a QCD constraint on only one of them. Such 
heavy quark symmetry relations are spoiled only by spin-symmetry violating corrections, 
which are expected to be smaller than the flavor-symmetry violating corrections inherent 
in the method of Sec. 3. 

All of the errors described above are either extremely small, or amenable to systematic 
reduction. We see no theoretical obstacle to predicting the B — ► D and B — > D* form 
factors to 1% accuracy, given the normalization at threshold and sufficiently precise mea- 
surements to fix two parameters. Such a prediction would not only test our understanding 
of QCD at an unprecedented level; it would present a precision probe of non-standard 
model physics in a hadronic arena. 



4-2. Reliability of the Extraction 

One type of heavy quark symmetry correction arises from the application of the B — > B 
dispersion relation to B —>■ D*lu. We estimate such corrections by making a 20% change 
in the ranges of ao and a±, resulting in a 2% shift in the central value of \V c b\. 

Another heavy quark correction arises at O(j^) in the normalization of the Isgur- 
Wise function at threshold. The normalization of the form factor g(uj = 1) has been 
estimated to be g(l) = 0.96JTJ, g(l) = 0.890, and g(l) = 0.93[|H§. We have included a 
QCD correction of 0.985, so to good approximation, this simply rescales the values of \V c b\ 
in Table 1 by 

There are other errors in our extraction that are not purely theoretical. The most 
pressing of these involve the binning of the measured rate against uj, smearing of uj in- 
troduced by boosting from the lab to the center of mass frame, and correlation of errors. 
Randomly varying input values of uj in our least squares fit of the CLEO data by ±0.05 
changes the central value of \V c b\ by less than 1%. A more thorough extraction can be 
done by the experimental groups themselves, using our basis function expansion in their 
maximum likelihood programs. 



5. Implications for \V u b\ 

The parametrization of form factors in terms of our basis functions applies to other 
heavy hadron decays, including B — > nW. In this case the range of the kinematic variable is 
larger, < z < 0.5, so more coefficients a n are needed for comparable accuracy. We expect 
six to eight a n will be necessary for accuracy of a few percent over the entire kinematic 
range, depending on the form of the actual data. 

A least squares fit will only produce best values of the products and |Kt&|a n , 

destroying hopes of a model-independent extraction of \V U b\- However, small values of \V u b\ 
tend to wash out the nontrivial z dependence, while the a n 's cannot compensate because 
they are bounded from above, so the extraction of a model-independent lower bound on 
\V u b\ should be possible. 

To obtain an upper bound on \V u b\ will require as input the overall normalization 
of the relevant form factor. Practically speaking, this means using a model or lattice 
simulation. Any candidate model must predict a form factor that is consistent with the 



basis functions, in the domain of validity of the model. This is a severe test to pass [15 
and should serve as an effective discriminator for models. 
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6. Conclusions 



The extraction of the CKM mixing parameter \V c b\ involves several types of uncer- 
tainties. Typically, these uncertainties are classified as 

\Vcb\ = V ± {stat} ± {syst} ± {life} ± {norm} ± {param} (6.1) 

where stat and syst refer to statistical and systematic experimental uncertainties, life refers 
to uncertainties in the B lifetime, norm refers to uncertainty in the value of the form 
factor at threshold, and param refers to uncertainty in the extrapolation of the measured 
differential rate to threshold. 

Not only the central value, but also the statistical uncertainty depends on the 
parametrization. For example, linear fits to CLEO data yield substantially smaller sta- 
tistical uncertainties than quadratic fits. Typically, quoted values correspond to the 
parametrization yielding the smallest statistical uncertainty, in effect throwing some sta- 
tistical uncertainty into the parametrization uncertainty, which remains implicit. Clearly, 
this does not improve the accuracy with which we know \V c b\. 

In this paper, we have essentially eliminated the uncertainty in the choice of 
parametrization. This was accomplished in four stages. First, we used QCD disper- 
sion relations to constrain the B — > B elastic form factor. Second, we derived a set 
of parametrized basis functions which automatically satisfies the dispersion relation con- 
straint, and expressed the B — > B form factor in terms of this basis. This expression 
involved an infinite number of parameters bounded by ^^L l°™| 2 — I- Third, we 
used heavy quark symmetry to relate the B — > B form factor to B — * D*W form factors 
and fixed the normalization at threshold. Over the entire kinematic range relevant to 
B — > D*W, we showed that neglecting all but the first two parameters ao, a± resulted in 
at most a 1% deviation in the predicted form factor. Finally, we made a least squares fit 
of the differential B — > D*W rate to |V c b|, ao, and a\. 

The results of this fit improve on all previous extractions in one important way: The 
uncertainty due to the choice of parametrization is under control, and of order 1%. Our 
statistical errors are larger than many quoted values. This does not reflect an inferiority 
of our method, but rather quantifies uncertainties that were previously left implicit. An 
averaged value from CLEO, ARGUS, and ALEPH data is 

\V cb \ = 0.037™ ( s tat). (6.2) 
9 



An estimation of systematic uncertainties requires a detailed knowledge of the experiments. 

The basis function parametrization described here allows generalization to B — > tvW 
as well. In this case, we expect to be able to extract a lower bound on \V u b\. In addition, 
precision tests of QCD-predicted form factors are now possible; these should be useful as 
checks of QCD models and lattice simulations. As experiments improve, they may even be 
used as probes of new physics. Applications and generalizations of the methods described 
here look promising. 
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